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A sufficiently massive collapsing star will end its life as a spacetime singularity. 
The nature of the Hawking radiation emitted during collapse depends critically 
on whether the star's boundary conditions are such as would lead to the eventual 
formation of a black hole or, alternatively, to the formation of a naked singularity 
This latter possibility is not excluded by the singularity theorems. We discuss 
the nature of the Hawking radiation emitted in each case. We justify the use 
of Bogoliubov transforms in the presence of a Cauchy horizon and show that if 
spacetime is assumed to terminate at the Cauchy horizon, the resulting spectrum 
is thermal, but with a temperature different from the Hawking temperature. 
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1. Introduction. 

There are by now many known examples of formation of naked singularities 
in spherical gravitational collapse in classical general relativity.' 1 ' 2 1 Whereas these 
examples do not necessarily invalidate the Cosmic Censorship Hypothesis, ® it 
is interesting to ask what a star forming a naked singularity would look like to 
a distant observer. Since visible regions of very high curvature develop during 
the collapse, it can be expected that quantum effects will play a significant role 
in determining the evolution of the star. Furthermore, it should be possible to 
describe these quantum effects using techniques from quantum field theory up to 
the time when curvatures approach Planck scales. 

A principal issue is a comparison between the Hawking evaporation of a star 
that forms a black hole, and the corresponding quantum evaporation of a star that 
forms a naked singularity. Studies of this problem can be divided into two classes; 
one in which the vacuum expectation value (VEV) of the energy momentum tensor 
of a quantized field is calculated in the background of the classical collapsing star 
using the trace anomaly, and the other in which the VEV of the radiation flux 
and spectrum of the radiation is calculated asymptotically in the geometric optics 
approximation, using Bogoliubov coefficients J 4 ' 

A few studies have been carried out in recent years to calculate the VEV 
of the quantized stress tensor in spacetimes which evolve to naked singularities. 
The central idea here is to investigate the behavior of the VEV in the approach 
to the Cauchy horizon. It has been found in all examples of shell focussing naked 
singularities that the flux of radiation diverges as the Cauchy horizon is approached. 
The divergence of the outgoing flux of the quantum field on the Cauchy horizon 
would suggest that the back-reaction ultimately prevents the naked singularity 
from forming. Furthermore, a divergent flux could in principle be measured by a 
distant observer if such objects were to occur in nature. 
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Perhaps the first investigation in this context was due to Ford and Parker, ^ 
who calculated the outgoing flux of a quantized massless scalar field in the space- 
time of a collapsing spherical dust cloud which develops a shell-crossing naked 
singularity. The flux was calculated using the geometric optics approximation, 
and remains finite in the approach to the formation of the naked singularity. An- 
other early calculation was by Hiscock et al.,® who computed the outgoing flux 
for a massless scalar field in the two dimensional self-similar Vaidya spacetime 
(obtained by suppressing angular coordinates in the spherical spacetime) which 
evolves to a shell-focusing naked singularity. In this instance, the flux diverges on 
the Cauchy horizon. In recent works' 7 ' we calculated the outgoing flux in the back- 
ground spacetime of a spherical self-similar dust cloud with a naked singularity, 
using both the geometric optics approximation (analogous to Ford and Parker' 5 ]) 
and the trace anomaly (analogous to Hiscock et al.' 6 '). Again, the flux diverges on 
the Cauchy horizon. 



As mentioned above, apart from the calculation of the quantum stress tensor, 
a calculation of the spectrum of the created particles is also of great interest. Since 
the calculation of the stress tensor is local, it can be carried out using standard 
techniques of quantum theory. However, the presence of the Cauchy horizon raises 
subtle issues in the calculation of the spectrum, and the purpose of the present 
paper is to address some of these issues. 



Consider the Penrose diagram for a collapsing star which develops a naked 
singularity (figure 1). 
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Figure 1. Formation of a naked singularity 

The presence of a Cauchy horizon implies that the region of X + to the future of the 
Cauchy horizon is exposed to the naked singularity. Hence one might conclude that 
a complete set of modes cannot be defined on Z + , and the standard Bogoliubov 
calculation cannot be carried out®. 

There are two ways out of this apparent obstacle. One is to consider the very 
real possibility that a star collapsing to a naked singularity does not destroy the 
whole universe. X + continues to be well-defined beyond the Cauchy horizon, and 
since no ingoing modes reach out to this part of I + , the outgoing modes in this 
region can be set to zero. It has been shown by us that the resulting spectrum 
is non-thermal.® In section 2 we give a justification for this approach by showing 
that the total radiation computed from this spectrum indeed equals the integrated 
flux obtained from the stress tensor. 

A second way out of the apparent obstacle is to assume that the Cauchy horizon 
is actually the end of spacetime, i.e., that the analytical continuation which defines 
the space time to its future is not physical. We show in section 4 that in this case 
the spectrum is black-body, with a temperature different from that for the Hawking 
black hole. This case is analogous to that of a marginally naked singularity (treated 
in ref. [6]), i.e., one in which the Cauchy horizon coincides with the event horizon. 

Another issue in the calculation of the spectrum concerns the existence of 
a complete orthonormal basis set of infalling waves on the Cauchy horizon. In 
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principle, the quantum field should be expressible in complete bases on X + U He, 
where He is the Cauchy horizon. There is an obvious difficulty in constructing a 
basis of infalling waves on He on account of the central singularity, and so there 
is an essential ambiguity in evolving the field in the future of the Cauchy horizon. 
This ambiguity is inconsequential, for it would be relevant if one were interested 
in constructing outgoing wave packets on X + in the future of the Cauchy horizon, 
from infalling waves on X~ . No such outgoing packets exist, however, because 
all probes coming in from Z~ at such late advanced times would be absorbed by 
the singularity. We are therefore only interested in complete orthonormal sets of 
solutions to the wave equation on hypersurfaces in the past of the Cauchy horizon. 
These are well defined. Moreover, the spectrum on X + will be independent of any 
particular choice of infalling basis states. 

We will compare the features of the Hawking radiation from black holes and 
naked singularities via a model of self-similar collapse based on inhomogeneous 
dust which is described in section 3. Our reason for using this particular model 
is twofold: firstly, the causal structure of the spacetime is well understood^ and 
secondly, the model exhibits both classical black hole and naked singularity end 
states, thus allowing for a comparison between the behaviors of each case within a 
unified picture. We do not expect any of the conclusions to be heavily dependent 
on the model itself, based as they are on general arguments that would apply 
whenever black holes or naked singularities are formed. 

2. Relationship between the flux and spectrum. 

Consider the propagation of a scalar field in some spherically symmetric, asymp-| 
totically flat background spacetime. For convenience we consider a massless scalar 
field, though massive fields as well as fields of arbitrary spin may be treated using 
the same techniques. We follow closely the treatment of Ford and Parker' 5 ! (see 
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also Birrell and Davies in ref. [4]). In this spacetime, we will assume that the 
radial null rays define a one-to-one mapping between a portion of past null infin- 
ity, X - , and some region of future null infinity, X + . We will also assume that a 
complete basis set of null infalling rays may be defined on X~ and a complete basis 
of null outgoing rays defined on X + . Let t,r,9 and define a quasi-Minkowskian 
coordinate system in the asymptotic region and let U = t — r and V = t + r be the 
the null coordinates there. We imagine that a null incoming ray, at V = constant, 
originates on X~ and propagates through the spacetime geometry, turning into a 
null outgoing ray, U = constant, on X + with value U = F(V). In the time reversed 
situation, one could trace a null ray, U = constant, on X + into the past. Such a ray 
would have originated at V = G(U) on X - , so that G(U) is the inverse of F(V). 
In Minkowski space, for example, F(V) = V and G(U) = U . If the functions, 
F(V), are known or can be determined, one considers positive energy solutions of 
the massless wave equation which have the following asymptotic form 

fulm - -^[e-^ + e-^lyUM), (2.1) 

which corresponds to an outgoing plane wave on X + and is normalized on a spatial 
hypersurface in the asymptotically flat "out" region according to 

(fulmifw'l'm') = 5(u - u')5u>5 m m', (2.2) 

with the inner product being defined by 

(7J) = -ij dZ"y/gZ [7(^7T) - {d~f)t] , (2.3) 

E 

E being the hypersurface. Wavepackets formed from the / w / m are outgoing plane 
waves at late times and incoming at early times in accordance with popagation by 
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geometrical optics, and any positive energy solution of the scalar wave equation 
that is outgoing on X + can be written as a wavepacket formed from the / w / m . An 
equivalent expression can be given for wavepackets that are incoming plane waves 
on X~ . These will be formed from 

Llm ~ -^[e-^ + e-^lyU^)- (2-4) 
A quantum field, 0, may therefore be expanded in either basis 

oo 

= 5^ J ' dujlaulmJulm + ^LlmTojlm] 
Im o 

oo 



(2.5) 



Im q 



in terms of the annihilation and creation operators, a^im and a u i m , and their 
hermitean conjugates. The vacuum defined by the a w / m , according to a^ZmlO) = 0, 
is the "in" vacuum and that defined by the a^/ m is the "out" vacuum. One is 
normally interested in the production of particles on T + , i.e., in the quantity 
(0|iV W 2 m |0) as V — > oo, where N^im — ^luflulm is the number operator in the 
"out" vacuum. It is then easily shown that this quantity is determined by the 
second of the two Bogoliubov coefficients, 

a LUU)' = (flVl fw') 

(2-6) 

that relate the two descriptions of the quantum field in (2.5), according to 



oo 

J dJ\(3{u'i 



\XM = J du'\(3{u'u)\ 2 . (2.7) 



We have suppressed the dependence on l,m because the geometric optics approx- 
imation that will be used in this paper is invalid for higher angular momentum 
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modes. However, these modes are expected to contribute little either to the flux 
or to the spectrum of the radiation because of the centrifugal potential, which 
would cause them to scatter to infinity before they encounter the region of high 
curvature(see, for example, B. S. DeWitt in ref.[4]) In the "out" region, 



f « 1 -iuG(U) 

7 ~ _ p -iuU 

Juj ~ r . e > 

V An w 



(2.8) 



giving 



U 



fttJu) = J dUe-^e-^\ (2.9) 

— 00 

where we have used (2.6), and where U = U represents the last outgoing ray 
that originated in an incoming packet from X~ . An equivalent and alternative 
expression, which constructs (5{uj'uj) on Z~, is 

= 1-^J dVe-^e-^ F{ y\ (2.10) 



using the asymptotic forms 



1 „-iwV 



1X0J 



1 « 1 -iuFjV) 
J u> ~ /-: e 



(2.11) 



on X , and where V" = V Q is the last incoming ray that turns into an outgoing 
packet on J + . 
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One may now compute the components of the stress energy tensor of the scalar 
field from the usual expression 

(0\T^(x)\0) = ljmV^lG^(x,x'), (2.12) 
x'— *x A 

where G^\x, x') is Hadamard's elementary function, 

G^(x,x') = (O|{0(x),#r')}|O>, (2-13) 

and 1)^' (x,x') is a non-local operator defined by the form of the stress-energy 
tensor and point-splitting. The expression in (2.12) must obviously be regularized 
and then renormalized, upon which one obtains, in particular, an expression for 
the radiated flux on X + J 4 ' 



F'" 3 / F" ' 1 
F^~ 2\F 



(2.14) 



24tt 

and a corresponding relation in terms of G{U). 

It is more interesting, however, to recover the radiated flux, (2.14), from a 
consistency condition: the integrated flux over I + must equal the total radiated 
energy as calculated by integrating the radiation spectrum over all frequencies, i.e., 

Uo oo oo oo 

J dU(0\T m \0) = J duu{0\N(w)\0) = J ' duuj J duj' (3*{uj'uj)(3{uj'uj). (2.15) 

-oo 

Using U = F(V), we can transform the integral over future null infinity to one 
over past null infinity, and write 

Vo OO 

J dVF\V)(0\T m \0) = J duu{0\N(u)\0) 

° (2.16) 

V V oo oo 

= JL J dV J dV> J 'dcv J ' du/u/e^-^eWW-^ . 

— oo — oo 
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While it may seem rather indirect, our reason for computing the radiated flux in this 
way is the following. When the singularity is globally naked, the Cauchy horizon 
will intersect future null infinity in the retarded past of the apparent horizon, at 
some point, say U Q . This means that the integration in expression (2.9) for j3{uj'uj) 
will not extend over all of X + , as it does for the black hole, but only the portion of 
X + that is in the retarded past of U — U . It is because, as we have pointed out in 
the introduction, any ray originating at such a value of V on X~ as would translate 
into an outgoing ray in the future of U — U would never arrive at X + , being, 
instead, absorbed by the singularity. If (2.16), for any V 0) reproduces the correct 
expression, (2.14), for the radiated flux as computed by a direct application of 
(2.12), it increases our confidence in the spectrum obtained from (2.7) even when 
the singularity is globally naked. V represents the last null ray that originates 
on X - and is able to reach X + . Both sides of equation (2.15) are infinite, but 
this does not trouble us as we are interested only in obtaining the flux. The 
infinite result is because, when the back reaction of spacetime is not accounted for, 
particle production will occur indefinitely even though the system loses energy. 
Energy conservation requires, therefore, that the back reaction will dominate at 
some stage. Before this stage is reached, however, (2.16) should serve as a good 
approximation of the actual physical situation. 

The integrals over oo and u/, on the right hand side of (2.16), can be performed 
to yield, 

V„ _ Vo Vo 

[ dVF'(V)(0\T Tlf M = I dV [ dV'^ = K (2.17) 

J 4tt 2 J J (V'-V) 2 \F(V')-F(V)} 1 ; 



— 00 — OO — 00 



The integral on the right has a pole at V' — V and contributions from points 
V' 7^ V vanish identically (as is seen by interchaging V and V' in (2.17)). Calling 
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z = V — V', the r.h.s. of (2.17) becomes 



Vo 



—I 

47T 2 



dV 



dz- 



V-Vo 



z 2 [F(V - z) - F(V)\ 



(2.18) 



We will define the z integral as the contribution from the portion of the contour, 
in the complex z— plane, shown in figure 2, that runs from V — V to — e along 
the real line (II), along the infinitesimal semi-circle, C , of radius e around the 
origin in the upper half plane, and from +e to infinity along the real line (III). The 
non- vanishing contribution from the z— integral in (2.18) is then just — irib-i(V) 
where where b-i(V) is the residue of the integrand at z — 0. 




Figure 2. The Contour C 

This is easily seen as follows: the value of the contour integral is identically zero, 
i.e., I + II + III + C = 0, as is the contribution from the semi-circle at infinity, 
Cqo. Thus we may write the integral in (2.18) as 



Vo oo 

dV I <!:- 



Vo 



V-Vo 



z 2 [F(V - z) - F(V)} 



dV / dz 



z 2 [F(V — z) — F(V)] 



(2.19) 



Next, consider the contour shown in figure 3. 
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Figure 3. The Contour C 

The value of this contour integral is 27ri&_i(V), i.e., I+II+III+C' = 2mb-i(y). 
Combining these two results, one finds I = itib-i(V) — (II+III). However, because 



V 



dV 



dz- 



-oo 



z 2 [F(V - z) - F{V)} 



(2.20) 



it follows immediately that value of the integral, as defined above, is 



Vo 



dV 



dz- 



1 



Vo 



4vr 2 J ' J " z 2 [F(V — z) — F(V)] 
-°° V-Vo 



-— / dVb-i{V). (2.21) 

4:71 



The residue can be evaluated by expanding the integrand about z = 0. We have 



z 2 [F(V - z) - F(V)} 



z F z z F- 



2! F' 3! F' (2!) 2 \F f 



+ 



giving 



b-i(V) 



6F' 



3 / F^ 
2\F> 



F 7 



(2.22) 



(2.23) 



or, inserting (2.23) into (2.21) and using (2.17), precisely the result in (2.14). The 
flux was originally obtained from (2.12) using standard point-splitting techniques. 
We have recovered it directly from a consistency condition, independently of details 
of the collapse. 
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3. The Collapse of Inhomogeneous Dust. 

We will apply the expressions obtained in (2.9) and (2.23) for the radiation 
spectrum and flux respectively, to the marginally bound, self-similar collapse of 
inhomogeneous dust. Although this model has been examined in detail by us 
elsewhere,! 6 ' 7 ' we include here a brief analysis of the causal structure of the space- 
time, both for the sake of completeness as well as to set our notation for the 
succeeding sections. It is described by the stress energy tensor 

7> = e(t,r)5X- (3-1) 

The metric is well known and given in comoving coordinates! 10 ' by 

ds 2 = dt 2 - R 2 (t,r)dr 2 - R 2 (t,r)dn 2 , (3.2) 

where the dust cloud is thought of as made up of concentric shells, each labeled by 
r. R'(t,r) is the derivative of R(t,r) with respect to r and R(t,r) is the physical 
radius (the area of a shell labelled by r is AixR 2 (t,r)) obeying, in the particular 
case of the marginally bound self similar collapse, 



R(t,r) 



3V\t 
2 r 



2/3 

(3.3) 



The physical radius is seen to depend on one parameter, A, (the "mass parameter"). 
This parameter determines the total mass, M(r), lying within the shell labeled by 
r as 2GM(r) = \r. The total mass of the dust is therefore 2GM = k = \r where 
r Q labels the outer boundary of the cloud. Now it can be shown that R(t, r) = is a 
curvature singularity. This means that the singularity curves are t (r) = 2r/(3VX), 
so that the last shell becomes singular at the time t a = 
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Beyond r = r Q spacetime is described by the the Schwarzschild solution 



ds 1 



(i-£)<«*-(i-£)"V-* w 



(3.4) 



and the first and second fundamental forms of the two patches must be matched 
at the boundary. This has been done in the past and gives 



To(t) 



- 2\ kR q - -R \ — + «Tn ^= — 

3 V « V^-v^ 



= t 



,3/2 



2y kR + /tin 



go + V^ | 

Ro- y/K 



(3.5) 



i? (t) 



-, 2/3 



1 — a- 



where we have set a = 3\/\/2. 

For the marginally bound, self-similar collapse under consideration, it is rela- 
tively simple to find null coordinates for this system. Consider the effective two 
dimensional metric, 



ds' 



dt 2 - R 2 (t,r)dr 2 , 



(3.6) 



and change variables to z,x where z = lnr, x — t/r. This gives 



ds 2 = r' 



dx 2 + 2xdxdz + (x 2 - R 2 {x))dz 2 



= r\x 2 -R%dr> - d X % 



i'2 



(3.7) 



where 



(3.8) 



in terms of 



I ± (x) 



dx 



x±R' 

We would like to choose null coordinates such that in the limit as A 



(3.9) 
these 



14 



reduce to the standard null coordinates in Minkowski space. Such coordinates are 
given by 



+re 



(3.10) 



x - R' > 
x - R' < 
x + R! > 
x + R' < 

To further analyze the causal structure, it is now convenient to go over to the 



u = 



V = 



-re 



+re* 



-re 



variable y defined by y = \ R/r. In terms of y, the integrals I± can be written as 



V dy 



3t/ 4 =F ay 3 - 3y T 2a 



(3.11) 



and the coordinates (3.10) become 



+re 



u 



— re 



7- 



where 



/-(y)<o 
/-(y)>o 

+re J + /+(y) < 
-re J + /+(y) > 



f±(y) = 3|/ 4 =F ay 3 - 3y T 2a. 



(3.12) 



(3.13) 



Let af be the roots of f±(y), for i e {1,2,3,4}. As f± are both real, they admit 



either 0, 2, or 4 real roots. The integrals can now be put in the form 



3 / dy 



E 



r (y - «f ) 



(3.14) 



where the A i are constants related to the coefficients of f±(y) and their roots by, 

.±3 



(3.15) 



1 ~ /i(af)' 

In particular, the Af satisfy J2i Af = 1- ^ an ^he roots are real, the solution is 
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explicitly given by 



u(y) 

v(y) 



± r Y[ \y - 



1 3 A; 



4 



(3.16) 



y 



a 



i=l 



We will now consider the case in which there are two real roots and a conjugate 
pair of complex roots. As we will shortly show at least two real roots (possibly 
degenerate) are required for the existence of a globally naked singularity at the 
origin so we do not consider the case when all the roots are complex even though 
it may be carried out in the same spirit. Let us order the roots so that the first 
two, ai ; 2, are a complex conjugate pair and 0:3,4 are real. From (3.15) it follows 
that A\ t 2 is also a complex pair whereas A^^ are real. Then the integrals are of 
the form 



I = 3 dy 



4 . 
Aj 

^ (y - «i) 



A\n(y-a) + A* \n(y - a*) + ^ Ai In \y - 

i=3,4 



a. 



(3.17) 

where a, a* are the complex roots and A, A* are the (complex) coefficients. Putting 



\A\e 1 *, 



y-a 



\y - «|e^, 



(3.18) 



so that the u, v coordinates have the explicit (and formal) solution 



u(y) = ±r\y-a-\ 6 ^ cos ^e- 6 ^^ s ' m ^U i ^ A \ y -a i 
v(y) = ±r|y-a+| 6 ^ + l co ^ + e -^ + ^ +si ^ + n^3,4b-a 



1 3 a; 



+ |3A+ 



(3.19) 



Consider the center (r = 0) at early times, t < 0. Then, because y — (1 



1G 



at/r) 1 / 3 — > oo, (3.19) gives (when all roots are real) 

|3E^r = _ r (i _ a l) ^ at 



u — > — r|y| 



* (3 ' 20) 

- r |y| 3 EiA + = _ r (l _ a _) _^ a t. 



r 



This line is therefore given by u — v. When two of the roots are complex conjugates 
of each other, the line is still u = v as we now show. Note that 



tan 1 



Im(— a) 



Re(y — a) 

(y is real), so that as y — > oo, £ — > 0. Then clearly 



u _, _ r | 2/ |3(2RcA-+A7+^) 



- r j/ 



3(2ReA++A++A+) 



(3.21) 



but, since ^ ^ = 1, we have the same result as before. 

The general solutions in (3.16) and (3.19) are useful to analyze another limit, 
namely the singularity at r — > at. This means that y — > 0. Now when y — > 0, 
/-(y) > and < 0. Then we see that (if all roots are real) 



u = - r | | \a { 



n 



(3.22) 



and, in particular, 



which is a negative constant, in general ^ — 1. The singularity is therefore spacelike 
until the last shell, r = r , collapses at £ = £ = r /a. The case of a pair of 
conjugate complex roots trivially gives the same result. Beyond this point the 
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singularity will be spacelike because it is just the Schwarzschild singularity in the 
exterior region. The behavior of the origin, r = 0,t — 0, is peculiar. It is the 
meeting point between two lines u = v and u = —cv and its nakedness (coveredness) 
is far from clear. However, if a null ray originating at this point reaches the 
boundary at Kruskal coordinate U < in the Schwarzschild region, it will reach 
X + and then the origin will be globally naked. 

We will be interested in the earliest null ray leaving the singularity and reaching 
1 + (the Cauchy Horizon) as well as the earliest null ray that strikes the singularity 
from X~ . These rays can be expected to intersect the first singular shell at r = 
0, t = 0, so it is natural to carefully examine the null rays passing through this 
point. The origin, being the intersection of the lines u = v and v = —cu (c/ 1 in 
general), corresponds to the point u — = v. Now any null ray traveling toward 
X + with u = must have either r = or I — > — oo. Therefore, when r ^ 0, such 
a ray is possible if and only if y = for some real root, of the polynomial 
f-(y). Indeed such a root may not exist, in which case the singularity is not naked 
as no null rays can emanate from it. In this black hole is formed, as shown 

in figure 3. 




Figure 4. Formation of a black hole 

If a real root exists however, at least one null ray leaves this point and reaches 
the boundary. The existence of real roots of the polynomials /_ (y) is therefore a 
necessary condition for the nakedness of the origin. This places a constraint on the 
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possible value of the constant a in the mass function. One finds that real roots exist 
provided that a < a c ~ 0.638 J 9 ' Each root corresponds to a null ray emanating 
from u = = v and there are at least two of them, if any at all. Because y = an 
implies that t = r(l — af)/a, we choose the largest real root of f-(y) as the one 
that gives the earliest null ray emanating from u = = v and call it cc_. Thus, 
y = a- is the Cauchy horizon. 

A similar reasoning can now be given for the incoming rays passing through 
u = = v . Again any ray with v = for r ^ must have i+ — > — oo, which 
is possible only if y = for some real root, of the polynomial f+(y)- Now, 
f+(y) admits two real roots, one unphysical (negative) and one positive. Again, 
call the (positive) physical root a+. 

What we have described above is pictured in the Penrose diagram of figure 1. 
It is to be expected that the behavior of quantum fields will be extremely sensitive 
to the collapse scenario being considered, that is to whether the mass parameter a 
lies below or above its critical value, and this is the topic of the next section. 



4. Radiation flux and spectrum for black holes and naked singularities. 

We will henceforth consider rays in the neighborhood of the lines given by 
y = a- for outgoing rays and y = a+ for incoming rays. The precise values of a± 
in terms of the mass parameter will not interest us for this work but we will Taylor 
expand about these two values, considering y± — y± + a±. 

Returning to (3.5), one can rewrite the Schwarzschild radial coordinate and 
time on the boundary as follows 

R (y) = r y 2 

T o{y) = -~y 3 - \ar y - \a 2 r In 1 ^^° \ \. 

a 3 9 3y/2a + 1 

Therefore, the Eddington-Finkelstein null coordinates on the boundary, U (y) = 
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T (y) — R *(y), V Q (y) = T {y) + R *(y), (where R * is the tortoise coordinate) take 
the form 

U (y) = - — y 3 - \ar Q y - r a y 2 - \a 2 r Q In \3y/2a - 1| 

4 8 ^ 

V {y) = -y 3 - -ar a y + r y 2 + -a 2 r In \3y/2a + 1|. 

do y 

It is now clear that the earliest null outgoing ray, u — 0, from the origin (the 

Cauchy Horizon) within the cloud strikes the boundary at y — a- and translates 

into the null outgoing ray 

= -a^_ — -ar a- — r a 2 _ — ^a 2 r In |3o;_ /2a — 1|, (4.3) 

which is never infinite (2a/ 3 is not a root of f-(y)). This null ray corresponds to a 
finite value of U and will therefore reach J + , so the existence of real roots of f-(y) 
turns out to be not just necessary, but a sufficient condition for the origin to be 
globally naked. The same argument applies to the infalling ray(s): the earliest null 
ray to pass through the origin is the ray corresponding to the value y — a+, or 

K (0) = -—a% - ^ar a + + r a+ + ^a 2 r In |3a+/2a + 1| (4.4) 
a 3 9 

and, again, since — 2a/3 is not a root of f+(y), V is not infinitely negative and 
such a ray will have come from J~. Thus, the existence of a positive real root of 
f+(y) is sufficient to ensure that at least one infalling ray from T~ will intersect 
the origin. 

The next question we must address is the relationship between the U, V coordi- 
nates in the exterior and the u, v coordinates (equations (18, 21)) on the boundary. 
This is difficult to do in general, but if we confine our study to rays that are "close" 
to u = and v = we can arrive at some conclusion regarding the quantum radia- 
tion on X + near the Cauchy horizon. "Close" will be taken to mean linearizations 
about y = a± respectively for incoming rays and outgoing rays. 
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First consider outgoing rays. For y ~ define y = y + a_ and find that for 
small y 

I ~ 7 _lny + C»(y), (4.5) 

where 

7- = (4.6) 

giving 

■u = — r|y| 7 — > y — a,- = y — J . (4.7) 
Therefore in terms of u (on the boundary) we can write U as follows 

U ~ C/(°)(a_) + r_(a_)(y-a_) = C/(°)(a_) + T_(a_) ^ , (4.8) 



where 



r a 3 



T_ = - 9— r < when a < a c . (4.9) 

a(3a- — 2a) 

Likewise, for incoming rays, put y = y + a+ and find that 

1+ = 1+ lny + 0(y), (4.10) 

where 

^ = TU^r (411) 

giving 

/ v \ 1/7+ 

v = -r|y| 7 + -> = (--J . (4.12) 

Thus, in terms of v (on the boundary) we can write V as follows 

V ~ V^(0)( a+ ) + r+(a+)(i/-a + ) = V^(a+) + F + (a+) ( --Y^ , (4.13) 



r, 
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where 



+ 



- 9 



r a+ 



a(3a+ + 2a) 



< when a < a f 



(4.14) 



We are now in a position to compute the radiated power close to the Cauchy 
horizon in the geometric optics approximation, consider a ray V = const, in 
the infinite past. We are interested only in the region on X + that is close to the 
Cauchy horizon, so the approximations in (4.8) and (4.13) will suffice. As the null 
ray crosses the boundary, we have 



v(v) = vW + r+l- 



V \ T+ 

r , 



(4.15) 



This expression can be inverted to give 



v(V) 



- r n 



V°-V 



7+ 



(4.16) 



where we have used the fact that T + is negative. Next, reflecting about the center 
(here, u — v) gives 

7+ 



U{V) = -T 



V°-V 



ir 



+i 



(4.17) 



Now as the outgoing ray crosses the outer boundary, we have the relation 



u(u) = - r_ 



u(v) = t/ (0) - |r. 



i 

U \ T- 

r o, 

v°-v 



(4.18) 



where now we have used the fact that T_ is negative. Thus, the right hand side of 
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(4.18) is F(V) and it has the form 



F(V) = A - B(V<® - V) 



~ 2+ 

7- 



(4.19) 



where B is a positive constant which is given in terms of the roots, a± given before. 
We can now write down the power radiated as a function of V using (2.14), 



TfjfjiV) 



A8nB 2 



l-Y 



7 2(y(0) _ y)2 7 



7^0, 



(4.20) 



where 



7 = ^ 

7- 



The result can be written in terms of U by inversion, 



(4.21) 



48tt 



l-Y 



Y(u(Q) - uf 



(4.22) 



which diverges as the Cauchy horizon is approached. 

Let us now consider the case when the origin is not naked, i.e., all roots of the 
polynomial f-(y) are complex. This means of course that the Cauchy horizon is 
formed in the retarded future of the event horizon as shown in figure 4. /,From the 
expression (4.3) for U, this implies that the event horizon intersects the boundary 
at y — > 2a/3. We will be interested in late times, so consider a ray that is close to 
the event horizon and that therefore intersects the boundary at y = 2a/3 +y. For 
such a ray, 

U ~ -4Mln|y|, (4.23) 

where M = 2arl/9 is the total mass of the cloud. Continuing backward into the 
cloud, it is necessary to retain only terms that are linear in y in the expression for 
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u. This ray translates into the ray 

2cz< 3 3 xj 
u = r exp(J_(— + y)) = r exp(-y) ~ r (H ), (4.24) 



or 

y ~ f(^-l), (4.25) 

which, when substituted back into the expression for £/ in (10), gives £7 as a function 
of u inside the cloud: 

U = -4Mln|^ ^-l) I- ( 4 - 26 ) 

Reflecting at the center (v = u), in terms of the advanced coordinate, v, within 
the cloud as 

U = -4Mln|^— - 1^ |. (4.27) 

We must now find a relationship between v and V in the exterior region. Let 
us suppose that the ray U — oo traced backward to the ray V = V . It is not 
important to know the precise value of V though this can be done and gives 
complicated expressions in terms of the roots of the polynomials f±(y). It is clear 
that this value, V of V corresponds to some given value y of y on the boundary. 
Consider a linearization about this value y (y = y + y) so that 

V = V + V'(y )y + ... (4.28) 

where V'(y) is the derivative of V w.r.t. y. Again, expanding v about this value, 
y of y on the boundary r = r Q gives 

v = r exp(I + (y + y)) = v + v'(y )y + ... (4.29) 

The precise values of V 0) v 0) V'(y ) and v'(y ) will not interest us for the following 
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analysis. What is important is that 

v = v + i^\(V-V ) (4.30) 
V'{y ) 

is linear to the order of interest, and that v = r (from (4.29)), so that 

U = -AM\n\ V ~ V ° \ = F(V), (4.31) 
B 

where B = Sr V / (y )/av'(y ) is an irrelevant constant. Applying (2.14), it follows 
that 

TfjfAU) « — s (4.32) 

uuK ' 192ttM 2 V ; 

to leading order. The radiation flux is seen to approach a constant as the horizon 
is approached. 

The marginally naked singularity does not exist in this model. The singularity 
would be marginally naked if the Cauchy horizon coincided with the event horizon 
for some value of a. However, we have seen that the event horizon is given by 
y = 2a/3, which is not a root of the polynomial f-(y) (except when a = 0). 
The singularity is therefore either naked (a < a c ) or covered (a > a c ), but never 
marginal. 

We now turn to the spectrum of the radiation emitted by the singularities. The 
famous black body radiation spectrum of the black hole (see Hawking in ref.[4]) is 
a direct consequence of the form of F(V), given in (4.31), and the fact that the 
integral in (2.9) extends over all of X + , so that all the outgoing basis states are 
sampled by the wavepackets formed from incoming plane waves that have scattered 
through the spacetime. 

When a naked singularity is formed, the scattering can occur arbitrarily close 
to the singularity, so that F(V) (as given by (4.19)) has a significantly different 
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dependence on the advanced coordinate V. But there is another and more im- 
portant difference. If we consider the possibility that a collapsing star does not 
eliminate all of the spacetime to its future, but that the spacetime continues as 
the analytic extension of the spacetime in the past of the Cauchy horizon and that 
J + continues to be well defined in its retarded future (and is therefore complete), 
then since no outgoing wavepackets formed from infalling plane waves are able to 
reach X + in the retarded future of the Cauchy horizon we see that not all of the 
outgoing basis states on J+ are sampled by the outgoing wavepackets. A direct 
consequence of this is that the spectrum is non-thermal. A simple calculation of 
the Bogoliubov coefficients, with F(V) given in (4.19), yields,® 



or 



The first expression, eq. (4.33) above, is useful to analyze the high frequency 
limit (uj'(Buj)^ 1 ^ — > 0) limit of the spectrum, for in this limit it is sufficient 
to consider only the first term in the series. Integration over uj' then yields the 
familiar logarithmic divergence and the spectrum is seen to fall of as 1/uj in the 
high frequency region. The second expression, eq. (4.34), serves to analyze its low 
frequency (u'(Buj)^ 1 ^ 1 — > oo) behavior. Integration over uj' in this limit shows a 
power law divergence in the infrared. This divergence is associated with the fact 
that there are an infinite number of quanta in each mode on X + . The difference 
between the divergence in the low and high frequency regimes may be associated 
with the red-shifting of modes in the proximity of the putative Cauchy horizon. 
Nevertheless, \{3(uj',uj)\ 2 is seen to be well behaved as a function of uj, falling as uj 
when uj — > 0. 
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As mentioned in the introduction, it is possible that the Cauchy horizon should 
be regarded as the natural end point of spacetime, so that the analytical contin- 
uation we have considered above is not physically acceptable. In this case, Z+ is 
not complete and U is no longer a good asymptotic coordinate. 

Referring back to figure 1, we see that the transformation to asymptotically 
flat coordinates must take the form 

U = -2Ke-°l 2K = -2Ke-^ 2K + , (4.35) 

which defines the asymptotic null coordinate u and where is defined by 

C/(0) = - 2 Ke-° m / 2K , (4.36) 

in terms of given earlier. Clearly, its definition is such that u ranges from 
— oo to +oo, while U ranges from — oo to the Cauchy horizon, A complete 

outgoing basis set will be defined w.r.t. u instead of U . Putting U = —2Ke u / 2n in 
(4.35) we find 

u = m - 2 K ln|l- e ^- f>( ° ) )/ 2K | ~ (/W - 2k1i. ° ' ' 



U - [/(O) = 2«;e-^ (0) )/ 2K 



2k 1 (4.37) 



(using the fact that we are near the Cauchy horizon). Now relating V and u 
(instead of U as we did earlier) we find 

u = F(V) = U - ^lnl^-^1, (4.38) 

7 D 

where B' is an irrelevant constant. However, this is precisely the relationship 
between the infalling and outgoing coordinates for a black hole, given in (4.31). It 
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will consequently yield thermal radiation, 

|/3(c/,cu)| 2 = —- A -. , (4.39) 

at the modified temperature, T = ^/Aim. This situation is analogous to the 
marginally naked singularity treated by Hiscock et. alJ 6 ] 

5. Discussion 

In this paper, we have used the marginally bound, spherically symmetric col- 
lapse of inhomogeneous, pressureless dust, which admits both classical black hole 
and naked singularity end states, to illustrate some key differences between the 
Hawking radiation from these objects. The central distinguishing feature appears 
to be the rapid flux of radiation that will be emitted from the naked singularity 
in the approach to the Cauchy horizon. The intensity of the radiation is clearly a 
consequence of the large curvatures that are encountered by infalling rays on their 
way out to future null infinity, so we expect this feature to hold true generically, 
and whenever regions of high curvature are visible to the asymptotic observer. In 
this spirit, the naked singularity may be thought of as a region of high curvature 
that is visible from future null infinity and not necessarily as a true singularity of 
spacetime. This rapid evaporation signals an instability of the Cauchy horizon, 
and may cause the collapsing star to evolve in such a way as to avoid its actual 
formation. It may be the mechanism by which nature avoids naked singularities, 
in which case the Cosmic Censorship hypothesis would originate in the quantum 
theory. We have not addressed the manner in which the appearance of the Cauchy 
horizon may be avoided as this requires a detailed study of the back reaction of 
spacetime. If nature does in fact employ the quantum theory to avoid naked sin- 
gularities, the magnified luminosity is likely to be observable and may allow for a 
glimpse into the behavior of matter fields in strongly curved backgrounds. Again, 
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if this possibility is taken seriously, it becomes necessary to look for additional 
features of the radiation that would distinguish naked singularities from other ra- 
diating objects. The spectrum of the radiation is one possibility. Although it is 
to be expected that the electromagnetic spectrum reaching the distant observer 
will not be characteristic of the collapsing star but rather of the thermalized debris 
surrounding it, the spectrum radiated in the form of neutrinos and gravitational 
waves should escape the surrounding matter relatively undisturbed. 

A central issue in the calculation of the spectrum of the Hawking radiation 
from spacetimes that admit Cauchy horizons, is that of the existence of a complete 
future null infinity. We know of no way to address this question within the semi- 
classical approach. There are, however, only two logical possibilities: (i) either the 
spacetime continues beyond the Cauchy horizon or (it) it terminates at the Cauchy 
horizon. 

If the local collapse of matter does not destroy the entire universe in its future 
(it is difficult to imagine that it would), the spacetime can be analytically continued 
beyond the Cauchy horizon and a complete future null infinity exists. Yet, no 
incoming waves on past null infinity are able to form wave packets to the future of 
the Cauchy horizon, therefore only a part of I + is actually probed. This results 
in a non-thermal spectrum. However, it also raises the issue of the consistency 
of the spectrum derived from the Bogoliubov coefficients. We have addressed this 
issue by showing explicitly that the total energy radiated, as computed from the 
integrated spectrum (derived from the Bogoliubov coefficients) is identical to the 
total energy radiated, as computed from the stress energy tensor to leading order. 
On the other hand, if the local collapse of a star indeed would destroy the universe 
in the future of the Cauchy horizon, so that the spacetime must terminate there, 
we showed that the spectrum of the radiation emitted is then necessarily thermal 
at a modified temperature. 



29 



Only a complete theory of quantum gravity can answer the question of the 
existence of a complete X + as this depends on the final fate of the collapse. We 
note, however, that the instability of the Cauchy horizon appears to signal that all 
of X + will survive and therefore that the spectrum will be non-thermal. 
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